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Perturbation theory can be reformulated as dynamical theory. Then a sequence of perturbative 
approximations is bijective to a trajectory of dynamical system with discrete time, called the approx- 
imation cascade. Here we concentrate our attention on the stability conditions permitting to control 
the convergence of approximation sequences. We show that several types of mapping multipliers 
O ' and Lyapunov exponents can be introduced and, respectively, several types of conditions controlling 

local stability can be formulated. The ideas are illustrated by calculating the energy levels of an 
anharmonic oscillator. 



m 



<D 



^ ' I. INTRODUCTION 



S 

-6 



Perturbation theory as applied to realistic physical problems practically always yields divergent sequences. There 
is a number of methods of finding effective limits of such sequences. The most known are the methods of Pade and of 
Borel summation. Constantly, new techniques appear. We are not aiming to give here a review of those. But we think 
it would be worth mentioning some of the recently developed interesting methods. Among them are some variants 
Q ■ of the Pade method [1,2] and the converse Cesaro method [3]. An analogy of the Kolmogorov perturbation theory 
O \ in classical mechanics has been constructed for self-adjoint operators [4], which reorganizes the asymptotic series in 
the usual Rayleigh-Schrodinger perturbation theory to a generalized asymptotic series. The accuracy of perturbation 
_^ | theory is known to improve when a perturbative expansion is accompanied by a variational procedure [5-7]. The 
£> . latter often gives good results already at the first step of calculations, as in the method of potential envelopes [8,9] 
0^ ' or when the zero approximation is very close to the considered problem [10]. The variational principle improves the 
accuracy of the quasiclassical and Langer methods [11] and of different sequence transformations [12]. The critical 
behaviour of several spin and vertex models has been studied by combining a variational series expansion and the 



o 



coherent anomaly method [13-17]. 

i- The common drawback of practically all perturbation techniques used for physical problems is the absence of a 

O"^ | general method for checking the convergence of perturbation sequences. This is because an explicit expression for 
the perturbative approximation of arbitrary order is, except a few trivial cases, never available. For some simple 
problems one can obtain exact numerical solutions. Then the accuracy of perturbation theory can be checked by a 
direct comparison of its results against the known exact answers. If, when increasing the approximation order, one 
gets better accuracy, then it is common to say that the procedure is convergent. Strictly speaking, this terminology 
is not correct, since the classical notion of convergence presupposes that the sequence of approximations has an exact 
solution as its limit when the approximation order tends to infinity. However, in actual calculations one always has a 
finite number of approximations, though the approximation order may be quite high. There is a number of examples 
• • ■ when the accuracy of a calculational scheme improves till a fixed approximation order but after it worsens resulting 
in a divergent sequence. Such behaviour is, indeed, a common feature of asymptotic series. In this situation it is more 
^^ , correct to say that the sequence of perturbation theory is semiconvergent. 

Note that here and in what follows we use the term " perturbation theory" in its most general meaning: perturbation 
theory is a regular procedure prescribing a general algorithm for defining a sequence of approximations of arbitrary 
order. The term "general algorithm" means that there is a general rule for calculating any approximation, although in 
practice it may happen that because of technical difficulties, one is able to find just a few initial approximations. The 
algorithm of perturbation theory may include, in addition to a basic expansion procedure, any of various resummation 
or renormalization tricks. In this sense, the so-called nonperturbative approaches are nothing but particular variants 
of perturbation theory, supplied with additional conditions. 

The main question for any variant of perturbation theory in the case of a complicated problem is how to control 
the convergence if neither an explicit form of high-order approximations nor the exact solution are available? To 
overcome this difficulty, an idea has been advanced [18] that the perturbation algorithm can be supplemented by a set 



of functions controlling the convergence of the approximation sequence. These functions because of their role, can be 
called the control, or governing, functions. Perturbation theory employing control functions has been first published 
in ref.[19] and used for describing anharmonic crystals [20-25]. In the cited papers, control functions were defined by 
a minimal-difference condition. A variational approach, called the minimal sensitivity condition [5], has also been 
applied to several anharmonic models [26-34]. The choice of conditions for control functions has been heuristic. 

To justify the option of conditions defining control functions, it was shown [35-39] that perturbation theory can be 
formulated as renormalization group theory. Then control functions are to be defined from a fixed-point condition, 
whose particular variants yield either the minimal-difference or minimal-sensitivity conditions. As far as a renormal- 
ization group can be considered as a kind of dynamical system, it was natural to reformulate perturbation theory to 
the language of dynamical theory [40-42] . This reformulation not only makes the theory more logical but also permits 
to define stability conditions related to the problem of convergence. In our previous papers we did not pay enough 
attention to these stability problems. The purpose of the present publication is to compensate this deficiency. We 
proffer a detailed analysis of the stability conditions showing that, really, there are several types of them, each with a 
different meaning. These conditions allow to understand intimate features of perturbation sequences and, therefore, 
to control their properties. Perturbation theory, whose convergence is controlled by control functions and which is 
supplemented by the stability conditions allowing for a thorough control of the properties of perturbation sequences, 
can be called the controlled perturbation theory. 

II. SURVEY OF APPROACH 

Before turning to stability conditions, we need to give a brief survey of the method whose stability is to be analysed. 

Assume that we are trying to solve a problem whose solution is a function f(g) , such that / : ED — > R , of 
a variable jgi C t . The first thing we have to do is to introduce control functions. For simplicity, we shall 
speak here about one set of such functions. The generalization to the case of several sets is straightforward [42]. 
Incorporate into the perturbation algorithm a parameter u £ R whose value is yet undefined. Then we get a 
sequence {Fk(g,u)}j^ of functions F& : ED x R — > R . In each term F k (g,u) , we substitute for the parameter u 
a function Uk{g) , such that Uk ■ ED — > R and k £ Z + = {0, 1,2,...} . This results in a function f k : D — ► R 
which is 

f k {g) = F k (g,u k (g)). (1) 

The role of the control functions u k {g) is to govern the convergence of the sequence {fk{g)} ( k Lo ^° ^ ne limit 

lim f k {g) = f{g). (2) 

k — >oo 

As is clear, there can be infinite number of different types of functions satisfying (2). This means that (2) defines a 
class of equivalence 

U = {u k {g)\k£Z+}. 

The latter can be explicitly found only if the general form of (1), as k — > oo , is known. Of course, this can be 
realised only for a few simple zero-and one-dimensional models [43-45]. Usually, the explicit form of (1) for arbitrary 
k — ► oo is not known. 

The second step is to narrow the class of equivalence U by those control functions that satisfy an evolution 
equation defining a dynamical system in discrete time, called the approximation cascade. This task is a kind of an 
inverse problem in the optimal control theory [46,47]. The direct problem in the latter is when an evolution equation is 
given and one has to find control functions minimizing a performance index. In our case, we need to find an evolution 
equation itself. To this end, introduce the coupling function g k (f) , such that g k : R — ► ED , given by the equation 

F (g k (f),u k (g k (/))) = /. (3) 

As is clear, the function g k x : ED — > R defined as 

g^(g) = F (g,u k (g)) (4) 

is inverse with respect to g k (f) , since 

9kH9k{f)) = f, 9k(g k l (g))=g. (5) 



Introduce the function 



Vk(f) = fk(9k(f)) (6) 



realizing an endomorphism y k : R — ► R of the measurable space R . The function (1) can be recovered from (6) 
by the transformation 

f k (g) = y k (9k 1 (g)). (7) 

The sequence {y/c(/)}^L ' by construction, is bijective to {fk(g)} ( j!° = o ■ 

Let us require that the endomorphism y k possess the semigroup properties 

Vk-y P = Uk+p, ya = 1- (8) 

This is equivalent to the relation 

Vk+ptf) = yk(y P (f)), (9) 

having the meaning of an evolution equation with the initial condition yo(f) = f . Equations analoguos to (9) can 
be met in various physical problems where they are often called the self-similar relations [48]. The requirement (8) 
narrows the class of equivalence IA to those control functions that provide the validity of the evolution equation (9). 
The semigroup 

Y = {y k } : Z+ x R -> R (10) 

of the endomorphisms y k defines a dynamical system with the discrete time k G Z + . In the dynamical theory this 
is called a semicascade. In our case, the latter is related to the sequence of approximations (7) because of which we 
call (10) the approximation cascade. The cascade trajectory {yk(f)}k=o * s bijective to the approximation sequence 
{/fc(ff)}/£o • The existence of the limit (2), which, according to (7), can be written as 

lim f k (g) = lim y k {g?{g)) - f(g), (11) 

k — *oo k— »oo 

is equivalent to the existence of an attracting fixed point of the approximation cascade, 

lim y k (f) = lim f k (g k (f)) = !/*(/). (12) 

k — >oo k — >oo 

At the third step, we embed the approximation cascade (10) into an approximation flow. This is done as follows. 
Instead of the discrete variable k G Z + consider a continuous variable t G IR+ = [0, oo) . Introduce an endomorphism 
y(t, •) : R — > R of the measurable space R , satisfying the semigroup properties (8) for each t G R+ . Then in the 
place of (9) we have 

y(t + t',f)=y(t,y(t',f)). (13) 

Require that the endomorphism y(t, •) would satisfy the conditions 

y(k,f) = y k (f) (fceZ+), 



The semigroup 



]im y(t, /) = !/*(/). (14) 

t—>oo 



Y={y(t,-)}: R+xR^R (15) 



of the endomorphisms y(t, •) is a dynamical system called semiflow. We name (15) the approximation flow. By 
construction, the semigroup (10) is a subgroup of the semigroup (15). In other words, the approximation cascade (10) 
is embedded into the approximation flow (15). The flow trajectory contains all the cascade trajectory, 

{yk(f)\keZ + }c{y(t,f)\teR + }. 



Assume that the embedding of the cascade into the flow is smooth, in the sense that the derivative dy/dt exists 
and is piecewise continuos on M + . Then the equation (13) can be rewritten in the differential form 



in which the velocity 



±y(t,f) = v(y(t,f)), (16) 



v(/) = lim^y(t,/) (17) 



is a transformation v : R — + T(IR) from R to a tangent space T(IR) . Integrating the evolution equation (16), we 
can cast it into the evolution integral 

y(t 2 J) d 

-f7=t 2 -*i. (18) 

The fixed point of the approximation flow is defined as a zero of the velocity 

«(»*(/)) =0. (19) 

The fourth step is the definition of quasifixed points, that is, of approximate fixed points. We cannot find fixed 
points exactly because the velocity of the approximation flow is not known, so we are not able to use (19). From 
another side, we cannot also define fixed points as the limit (12), since an expression for yk(f) at arbitrary k — > oo 
is not available. 

To find approximate fixed points, or quasifixed points, we need to have an explicit form of the cascade velocity 

v k (f) = v(y(t, f)) (t€[k,k + l}). (20) 

We may define it as an Eulcr discretization of the derivative dy/dt . In doing this, we take into account that Fk(g, Uk) 
depends on k directly as well as through Uk ■ In this way, we may write the cascade velocity (20) as 

v k (f) = V k (g k (f),u k (g k (f))) (21) 

with the finite difference 

d 
V k (g,u k ) = F k+ i(g,u k ) - F k (g,u k ) + (u k +i - u k )-^ — F k (g,u k ) 

ou k 

and the coupling function gk given by the constraint (3). As the cascade trajectory approaches the fixed point, when 
k — ► oo , then v k (f) — > . But to check whether the cascade velocity tends to zero, as k — > oo , we need to know the 
form of F k {g, u k ) for arbitrary k — > oo . If the latter would be known, then we could check directly the convergence 
condition 

Fk +p (g,uk+p) - Fk(g,u k ) — ► (k — > oo), (22) 

in which p > 1 . This condition would define those control functions which provide the convergence for the sequence 
{fk{g)} < k=o of terms (1). Such a condition of defining control functions may be called the asymptotic fitting condition. 
It has been used for some anharmonic models [43-45]. However, the possibility of finding the general expression for 
Fk(g,Uk) at the arbitrary k — > oo is rather an extreme exception. Usually, just a few first terms are available only. 
Therefore, we need to have some working conditions for defining control functions in the case of finite numbers k of 
F k . 

One such condition for defining control functions is the minimal-difference condition [18-25], 

F k+ i(g,u k ) -F k (g,u k ) = 0. (23) 

This condition makes zero only a part of the cascade velocity (21). Therefore (23) may be named the quasifixed point 
condition. 

Another quasifixed-point condition is the minimal sensitivity condition [26-34] , 

^-F k (g,u k )=0, (24) 

OUk 



which also makes zero only a part of (21). 

A slight generalization of (24) following from (21) is the quasifixed-point condition [42] 

d 

(u k+1 -u k )- — F k (g,u k ) = 0. (25) 

ou k 

The meaning of (25) is as follows: if (24) has a solution for u k , then this solution gives the control function u k (g) ; 
when (24) has no solution for u k , then we put u k = u k +\ . 

Defining the control functions from one of the quasifixed-point conditions, we obtain a sequence of quasifixed points 
Uk(f) and, respectively, a sequence of their images f k (g) that are the sought approximations for f(g) . Several 
other types of quasifixed-point conditions have also been studied [49,50]. 

In the fifth step of the considered perturbation theory we find corrected approximations. As far as the quasifixed 
point conditions do not make the cascade velocity exactly zero, the trajectory does not stop at y k (f) , though the 
motion slows down. If we accept the quasifixed-point condition (25), then the cascade velocity (21) in the vicinity of 
a quasifixed point y k (f) becomes 

vt(f) = V k *(g k (f),u k (g k (f))), (26) 

where 

V k *{g,u k ) = F k+ i{g,u k ) - F k (g,u k ). 

The motion in the interval of t G [k, k + 1] is described by the approximation flow (15). Substituting into the 
evolution integral (18) the time limits t\ = k and ti = k + 1 , we have 

Vk+iU) j,, 

Using in (27), instead of v(y) , the velocity (26), we get the evolution integral 

rvt(f) dy 



Vk(f) u k 



K(y) 



= 1 (28) 



defining the corrected quasifixed point y k (f) ■ Making in (28) the substitution / — > g k (g) , we come, according to 
(7), to the integral 

rfiia) df , N 

(29) 
J Ma) u k\J ) 

in which 



is the corrected k -order approximation. 



J Ms) v *k(f) 
r k {9)=yt{9k\9)) (30) 

III. STABILITY OF CASCADE 



A very important feature of the controlled perturbation theory is the possibility to control whether we are approach- 
ing the fixed point, that is, the correct answer, even if this exact answer is not known. This possibility is based on the 
semigroup property (9) of the approximation cascade, according to which each point y k (f) is considered as a result of 
mapping from a previous point y P (f) with p < k — 1 . We are approaching an attracting fixed point if the mapping 
is contracting. This, of course, depends on the perturbation algorithm and on the initial approximation. A given 
perturbation algorithm has a basin of attraction. The mapping can be contracting only if an initial approximation is 
in the basin of attraction. In iteration theory, the analog of the basin of attraction is the set of normality or Fatou 
set, while the complement of the Fatou set is called the Julia set. Attracting fixed points are in the Fatou set, while 
repelling fixed points are in the Julia set. When the iteration is done by means of an entire transcendental function, 
then every point in the Julia set is a limit point of repelling periodic points, that is, the Julia set is the closure of the 
set of repelling periodic points [51]. 



Let us analyse when the mapping corresponding to the approximation cascade is contracting. Consider the change 

y P (f) ^ y P (f) + Sy P (f) (31) 

with 8y p (f) — ► and the variation 

5y kp (f) = y k (y p (f) + Sy p (f)) - y k (y p (f)). (32) 

A particular case of (32) is 

Sy k (/) =y k (f + Sf) -y k (f) = 5y kQ (/) , (33) 

when the initial condition is changed. 

To proceed further, we need to introduce the notation for mapping multipliers. Define the quasilocal multipliers 

^(/) = -^jjy- ( 34 ) 

and 

Mfe (/) = M^ =/Ufe0 (/). (35) 

They satisfy the relation 

/%>(/) = Mfe(j/ P (/)) (36) 

and have the property 

/**(/) =/*(/) = 1. (37) 

Other useful relations can be derived basing on the semigroup property 

yk+p(f) = Vk(y P (f)) = y P {yk(f)) 



and the variational derivative 



In this way we obtain 



and 



Introduce the local multiplier 



fyk(f) = dy k (f)/df _ n k (f) 
$y P (.f) dy p (f)/df M P (/)' 

M/) - ^y^ (38) 

/%>(/)M/) = Mpfe(/)Mfc(/)- (39) 



^ (/) "^S)^ lfc - l(/) - (40) 



Using (38), this can be also written as 



/**(/) = -^TTT (* > !)• (41) 



In the case of fc = 1 , 

/»!(/)= Mi (/)• 

G 



The quasilocal multipliers (34) and (35) can be presented as products 

k+p 

tikp(f) = n ^u), 
j=p+i 

k 

»k(f) = Hn* p (f) (42) 

p=l 

of the local multipliers (40). In another form (42) reads 

l*kp(f) =Mfe+ p (/)Mfe-ip(/), 

Hk(f) = uUf)Hk-i(f). (43) 

With these multipiers, the variation (32) becomes 

Sykp(f) = Ukp(f)6y P (f). (44) 

Eq.(44) describes the deviation of the cascade trajectory at the (k + p) -step resulting from the variation Sy p (f) at 
a p -step. The mapping, corresponding to (44), is contracting if the condition of quasilocal contraction 

!/%>(/)! < 1 (45) 

holds. This shows that the mapping is effectively contracting after k steps starting from a p -step. Equivalcntly, 
one may say that the mapping is effectively contracting on the interval [p, k + p] . 

In particular, for the interval [0, k] , we need to deal with the variation (33) which yields 

Sy k (f) = n k (f)6f. (46) 

The condition of quasilocal contraction on the interval [0, k] is 

M/)l < 1, (47) 

where (J>k{f) is the quasilocal multiplier (35) with k > 1 . 

If we are interested in the contraction property at just one step, from k — 1 to k , then we have to consider the 
variation 

SVk(f) = l*Uf)SVk-i(f), (48) 

in which H%{f) is the local multiplier (40). The mapping is locally contracting at a A: -step if the condition of local 
contraction 

l/4(/)l < 1 (49) 

is valid. 

The condition (49) is stronger than (47) in the following sense: if \/jl\ < 1 holds for all p <G [0, k] , then (47) follows 
from this because of the relation (42), though the inverse is not true. It may happen that (47) holds, but for some 
p from the interval [0, k] the condition (49) is not valid. In other words, there can exist the effective contraction on 
an interval [0, k] although not for all steps there can be the local contraction. Simbolically, the relation between the 
two notions is: 

local contraction — > quasilocal contraction. 

The contraction for a mapping is the same as the stability for a cascade. The stability is characterized by Lyapunov 
exponents. Again, different kinds of such exponents can be defined. The quasilocal Lyapunov exponent 

A fep (/) = iln| Mfep (/)| (50) 

is related to the quasilocal multiplier (34), and the quasilocal exponent 



A fe (/) = -!- In M/) I = A fe0 (/), (51) 

to the multiplier (35). Taking into account (38), we have 

Xkpif) = ^A fe+P (/) - |A P (/). (52) 

The local Lyapunov exponent 

A£(/)s In |^(/)|=Ai *_!(/) (53) 

is defined through the local multiplier (40). Because of (42), the quasilocal exponent (51) is an arithmetic average 

fc 
k 



A t (/) = TA;(/) (54) 



P =i 

of the local exponents from (53). 

The effective stability on the interval [p, p + k] with p > and fc > 1 , means that the condition of quasilocal 
stability 

Xk P (f) < (55) 

holds. In the case of p = , this reduces to 

A fe (/)<0. (56) 

As is evident, (55) and (56) follow from the contraction conditions (45) and (47), respectively. 
The condition of local stability at a fc -step, from fc — 1 to fc , reads 

A£(/)<0, (57) 

resulting from (49). Anew, the condition (57) is stronger than (56) in the sense that the stability can exist on an 
interval but not necessarily at all points of the latter, while if (57) holds for all points of an interval, then (56) follows 
for this interval. 

The maximal Lyapunov exponent, usually employed in dynamical theory, is 

A(/) = lim Afc(/). (58) 

k— »oo 

The condition of asymptotic stability implies that 

A(/) < 0. (59) 

If the condition (56) is valid for all fc > 1 , then it is stronger than (59). Thus, the relation between the different 
types of stability is as follows: 

local stability — > quasilocal stability — > asymptotic stability. 

Recollect that the approximation-cascade trajectory {yk{f)}kLo is, by construction, bijective to the approximation 
sequence {fk(g)}kLo ■ Each point yt{f) has its image /&(<?) given by the relations (6) and (7). For the mapping 
multipliers and Lyapunov exponents introduced above, we may also define their images as functions of g . 

The image of the quasilocal multiplier (35) is 

M k (g)=u k (g^(g)). (60) 

This can also be written as 

»r / s flP ( , » SF k (g,u k {g)) 

°Fo(g, uk{g)) 

The image of the local multiplier (40) is 



M* k {g) = l4(9k-M)- (61) 

Using the properties of multipliers, we may write 

,.*/ n (r i w ^fc(-Fo(ff,Mfc-i(g)) 

M k (9) =Mi (/fc-i 3 = 7F7 ttt- 

The contraction conditions (47) and (49) can be reformulated for the multipliers (60) and (61) giving 

\M k (g)\ < 1, \M* k {g)\ < 1, (62) 

respectively. 

For the image of the quasilocal Lyapunov exponent (51), we get 

A k (g) = X k {gl\g)) = \\n \M k (g)\, (63) 

and for the image of the local Lyapunov exponent (53), 

K(g) = K(9 k li(9)) = ^ \M* k {g)l (64) 

These are connected with each other, equivalently to (54), through the arithmetic averaging 

k 

k 



k 

Mg) = lj2K(g)- 



The stability conditions (56) and (57) can be written for the Lyapunov exponents (63) and (64), so that 

A fc ( ff ) < 0, A* ( 5 ) < 0, (65) 

respectively. 

As far as the stability conditions for the Lyapunov exponents are a reformulation of the contraction conditions for 
the mapping multipliers, in what follows we shall often refer to any of them as to the stability conditions. 

IV. STABILITY OF FLOW 

The stability of motion analized above concerns the stability of an approximation cascade. If the latter is embedded 
into an approximation flow, we need to check the stability of the flow as well. This can be done by making in the 
evolution equation (16) 

the substitution 

y(t,f)->y(t,f) + Sy(t,f) (66) 

implying that Sy(t, /) — > . Then we find 

Sy(t, f) = Sy(t , /) exp U X(t', f)dt'\ , (67) 

where to < t , and 

is the local Lyapunov exponent for the flow. According to (20), the flow velocity in the time interval k < t < k + 1 
is given by the corresponding cascade velocity. For this reason, we put to — k in (67) and get 

6y(t, J) - %(/) exp{A fc (/)(* - k)}, (69) 

where we took into account that, in compliance with (14), y{k, /) = y k {f) , and that (68) gives 



A, (/) - %%P- (70) 

oyk(f) 

The approximation flow near a quasifixed point yk{f) is stable if 

Afe (/) < 0. (71) 

The image of (70) is 

Afe (5) =A fe (g k (3)) = s -, -t , (72) 

and the condition (71) of the local stability of the flow acquires the form 

Ak (g) < 0. (73) 

Note that the local Lyapunov exponents for a cascade and for a flow, into which the cascade is embedded, are, 
generally, different. This means that (53) does not coincide with (70). Therefore, it may happen that the cascade 
at a point yu{f) is locally stable but the flow at the same point is not, or vice versa. To understand better the 
distinction between (70) and (53), we may invoke the definition of the cascade velocity as of the finite difference 

Vk(f)=Vk+i(f)-Vk(f). (74) 

Then (70) becomes 

Afe (/) = Mfe+i fc(/) - Mfefe(/)- (75) 

Using the properties of the quasilocal multipliers (34), we may transform (75) to 

Afe (/) = nkkif) [iA fc+iC/O - 1] • 

According to (50), (51) and (53), we have 

{2fe 
p=fe+l 

| M fe(/)|=exp{fcAfe(/)}, 

K(/)|=cxp{A£(/)}. 
As we see, there is no simple relation between (75) and (53). 

V. STABILITY AND CONVERGENCE 

Since the cascade trajectory is bijective to the approximation sequence, the stability conditions for the cascade 
should characterize the corresponding convergence properties for the sequence. 
The deviation of the trajectory point yk(f) from the fixed point y*(f) is 

Ay k (f)^y k (f)-y*(f). (76) 

Consider the deviation Ayj. +p (/) assuming that y p (f) is close to y*(f) in the sense that 

lw P (/)-V*(/)l«IV*(/)l- (77) 

Employing the definition of the fixed point, 

y k (y*(f)) = y*(f), 
10 



we find 

Ay k+P (f) ~ n kp (f)Ay p (f). (78) 

From (78), it is easy to derive the equations 

Ay k (f) ~ M fe(/)A/, 

&y k (f) * l4(f)*Vk-i(f), (79) 

in which A/ = / — y*(/) . In its turn, (79) gives 

|Aife(/)|~|A/|exp{fcA fc (/)}, 

|A» fc (/)|~|A» fc _ 1 (/)|exp{AJ(/)}. (80) 



The image of (76) is 
Respectively, (80) yields 



Af k (g) = f k (g) - f(g). (81) 

|A/ fe ( 5 )|~|A/ (s)|cxp{fcA fc ( 5 )}. 

|A/ fc (s)| ~ |A/ fc _!( ff )| exp{A£(<?)}. (82) 

The accuracy of an approximation fk(g) , as compared to the exact value f(g) , is characterized by the absolute 
error (81). When the condition of quasilocal stability (56) holds, that is when A. k (g) < , then (82) shows that 
the accuracy of f k {g) is better than that of fo(g) , although the accuracy may decrease with respect to f k _\(g) . 
If Afe(g) < for all k > 1 , then the error (81) tends to zero, which means the convergence of the sequence 
{/fc(<?)}feLo ' When the conditions of local stability (57) is valid, then A£(.g) < , and (82) shows that the accuracy 
°f fk(g) is better than that of fk-i(g) ■ If the condition of local stability occurs for all fc > 1 , then the accuracy 
of approximations improves at each step. This also means the convergence of the sequence {fk(g)}'kLa > since (57) is 
stronger than (56). 

Thus, the stability conditions (65) are sufficient conditions for the convergence of {fk(g)}kLo ■ The necessary and 
sufficient condition of convergence is 

kA k (g) -> -co (fc -> oo). (83) 

For the corresponding quasilocal exponent (51), this reads 

fcA fe (/)^-oo (fc-oo), (84) 

which is equivalent to the asymptotic condition 

|Wt(/)|-0 (fc-oo). (85) 

The convergence conditions (83)-(85) are weaker than the conditions of quasilocal stability (56), local stability (57) 
and asymptotic stability (59). The maximal Lyapunov exponent (58) can be zero; nevertheless, the convergence of 
the approximation sequence will persist provided (84) is valid. 

By analogy with the usage of the terms quasilocal or local contraction, as applied to a mapping, and quasilocal or 
local stability, as applied to a cascade, we may use the terms quasilocal or local convergence, as applied to a sequence. 
We shall say that a sequence of approximations is quasilocally convergent on the interval [0, fc] if the condition of 
quasilocal stability (56), i.e., A k (g) < , holds on this interval. The sequence will be named locally convergent at a 
fc -step if the condition of local stability (57), that is, A* k (g) < , is valid at this given fc . When the conditions of 
quasilocal and local convergence hold everywhere, that is, are true for all fc > 1 , then they are stronger than the 
convergence criterion (63). There is the following relation between different notions of convergence: 

local convergence — * quasilocal convergence —* convergence. 

This notions should not be confused with the point convergence, convergence on an interval and uniform convergence 
of the sequence {fk(g)}kLo with respect to the variable g. 
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Because of the exponential renormalization of the deviation in (82), one may say, when the corresponding conver- 
gence conditions are valid, that there occurs an exponential convergence. 

In order to find out when the corrected approximation f k (g) is better that fk(g) , consider the deviation 

Ay* k (f) = y* k (f) -y*(f). (86) 

From the evolution equation (16) we obtain 

A^(/)~Ay fe (/)exp{A fe (/)}, (87) 

where Afc (/) is defined in (70). The image of (86) is 

A/ fc *(ff) = f k {g) - f{g). (88) 

Whence, (87) gives 

\Af* k (g)\~\Af k (g)\eMAk(g)}- (89) 

Eq.(89) shows that the corrected approximation f k (g) is more accurate than fk(g) when the approximation flow is 
locally stable, so that condition (73) holds. 

Let us also observe when the corrected approximation f k (g) is better than fk+i(g) ■ This happens when 

\Ay* k (f)\ < \Ay k+1 (f)\. 

The latter inequality leads to 

A* (/) < AJ +1 (/), 

Ak (g) < K +1 (g). (90) 

When the approximation cascade is locally stable at the (fc + 1) -step, then the improvement of the accuracy for 
f k (g) , as compared to fk(g) , can be achieved only if the enveloping approximation flow is stable at its k -point 
with the local exponent satisfying (90). If the approximation cascade is locally unstable, so that A£ +1 (/) > , then 
the corrected approximation f k (g) can be of much better accuracy than fk(g) even in the case of an unstable 
approximation flow, provided (90) is valid. In the latter case, to easier satisfy (90), the motion should be damped, so 

that to make Afc (/) smaller. This can be done by incorporating into the definition of the cascade velocity (21) a 
damping parameter 8k lowering Vk{f) , 

Vk(f) -> faVk(f)- 

The value of the damping parameter 8k can be found from additional conditions. For example, one may require the 
coincidence of some asymptotic values for the corrected approximation f k (g) and for the exact f(g) , of course, 
if such asymptotic values of f(g) are available [39,49]. Another option is to put 8k = \ , which corresponds to 
diminishing twice the step of the calculational procedure. 

The possibility of improving the accuracy even for unstable approximation cascades and flows is the main advantage 
of the corrected approximations. 

VI. ANHARMONIC OSCILLATOR 

To illustrate the ideas of the approach we choose an anharmonic-oscillator model. The anharmonicity of oscillations 
plays a very important role in many physical problems, for instance, in anharmonic crystals [52,53]. 
Suppose we need to find the energy levels of an anharmonic oscillator with the Hamiltonian 

in which x G (— oo, oo) and the coupling, or anharmonicity, parameter g > . Note that several problems of 
quantum mechanics can be reduced to oscillator-type models by a special change of variables [54,55]. 
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Emphasize that our aim here is not simply the calculation of the energy levels but the demonstration how the 
controlled perturbation theory, formulated as dynamical theory, works for such a touchstone model as (91). In our 
previous papers [38-40] we have considered solely the first step of the theory as applied to an anharmonic oscillator. 
One step, of course, does not permit yet to illustrate and exploit in full the underlying ideas. This is why we have 
again to turn our attention to the model (91) extending the consideration to the higher-order approximations. 

It is natural to start from the harmonic oscillator whose Hamiltonian 

contains an unknown parameter u . For convenience, we introduce the notation 

E k (g,u)=(n+^\F k (g,u) (93) 

for the k -order approximation of the spectrum. The quantum index n = 0,1,2,... in E k and F k is not written 
explicitly for the sake of brevity. 

The sequence {-Ffc(5,u)}^L } is to be obtained by the Rayleigh-Schrodinger perturbation theory starting from 



In what follows we shall need the notation 



F (g,u)=u. (94) 

a = a(u) = 1 -, 

(3 = 0(u) = ^f, (95) 

_ n 2 +n+ 1/2 
7= ra+1/2 

The first four approximations, under a fixed u , are 

u 
F 1 (g,u) = u--{2a-/3), 

F 2 (g, u) = Fiig, u) - | (a 2 - 2a/3 + 2af3 2 ) , 

F 3 (g,u) - F 2 (g, u) - ^ (a 3 - Aa 2 + I0aaf3 2 - 36/3 3 ) , 
F 4 (g, u) = F 3 (g, u) - ^- ( |a 4 - 8a 3 d + 35aa 2 d 2 - 24ba[3 3 + 4c/? 4 J , (96) 



in which 



b = 6( 7 ) 

c = c( 7 ) 



. , 17n 2 + 17n + 21 
a = a(7) = Wy) 2 ' 

_ 125n 4 + 250n 3 + 472n 2 + 347n + 111 
= (n + l/2)(6 7 ) 3 ' 

10689n 4 + 21378n 3 + 60616n 2 + 49927n + 30885 



8(6 7 ) 4 
From the quasifixed-point condition (25) of the form 

— F 1 (g,u 1 ) = 

OUi 
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we get the equation 

ul - mi - 67.9 = (97) 

for the control function u\(g) . Eqs.(95) and (97) tell us that 

a{ui) = P(ui). 
Introducing the notation 

a k = a(wfe) = 1 o ' Pk = /3(«fe), (98) 

we have 

F fc+1 (.g, Ul )= J F 1 fc ( 5 , Ul )+ Ul A lfe ^ +1 (As = 0,1), (99) 

where 

^10 = -|, A n = i (1 - 2a) . 
The equation dF2(g,U2)/du2 = has no real solutions for ui , therefore, according to (25), we put 

U2(g) = u 3 (g), (100) 

and u 3 being defined by the equation 

^3(^3) = 0. 

The latter yields 

u 3 3 - u 3 - 673.9 = (101) 

with 

73 = K7, a 3 = Kf3 3 (102) 

and with k given by the equation 

5k 3 - 24k 2 + 70a« - 246 = 0. (103) 

Eqs.(101)-(103) make it possible to find 

F k+1 (g,u 3 )=F k (g,u 3 )+u 3 A 3k a k 3 +1 (k < 3), (104) 



where 



.30 = 4 2-1 

4 V k 



A 31 = -\(l-l + 2a 3 ) 



1 / 4 
^32 = --7 l-- + 10a 3 -36 3 

16 V « 



^33 = -32 ( 4 " - + 35a 3 - 24& 3 + 4c 3 I . 
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and 

«3 = a(73), ^3 = b(l3), c 3 =c(j 3 ). 

Definition (3) of the coupling function gk(f) because of (94), gives the equation 

Uk(9k(f)) = }■ 
From here we find 

f(f 2 - 1) 



.%(/) = 



67/c 



where we took into account that 



and used the notation 



For the cascade velocity (26) we get 



with 



The evolution integral (29) becomes 



a{u k (g k (f))) = 1-72 



71 = 7> 72 = 73 = n-y. 



fe+i 



v * k (f) = A k f[l-- 



A k = A kk , A 2k = A 3k . 



f£(g) f2fc+l 



f k+1 df 



= A k . 



With the notation 



integral (109) transforms to 



I Mo) (/ 2 -l) fe+1 
ft (g) = yjl + z* k (g), f k (g) ee y/l + z k (g), 



z t(s) 



Zk(g) 



(l + zf 
yk +i 



dz = 2A k 



Employing the binomial formula 



(l + z) fe = ]TC^; C* S 

p=0 



it! 



(k — p)\p\ ' 



we integrate (111) obtaining the equation 



'fc-i 



z k = z k 



CI' 



ex p Err 



,p=0 



k — p 



1 1 



{z * )k - P z k- P 



2A k }, 



in which 



4=4(9), z k = z k {g). 



(105) 



(106) 



(107) 



(108) 



(109) 



(110) 



(111) 



(112) 
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k— p 

k — p 



Introducing the polynomial 

fc— l /~i%i 

with 

Pi(x) = x 



P 2 (x) = 2x+-x 2 , 



P 3 (x) = 3x + -x 2 + -x s , 
we can cast (112) into a more compact form 

zt = z k c^[p k f^J-P k f^-J+2A k \. (113) 

In this way, the k -order approximation for the spectrum of the Hamiltonian (91), i.e. 

e k (g) = E k (g,u k (g)), 
owing to notation (93) can be written as 

e k {g) = (n+-jf k (g) = e k {n,g) (114) 

with f k (g) = Fkig^Ukig)) ■ The corrected k -order approximation is 

«(.)-(» + £)/«,) -«(«..). (115) 

The accuracy of these approximations is characterized by their percentage errors 



£k{g) = r~, 100%, 

e(g) 
e%{g) = et{9) - e{9) • 100%, 

e{g) 

with respect to exact numerical values e(g) . 

Table I illustrates the accuracy of the approximations 

ei(g) = E 1 {g,u 1 {g))\ e 2 (g) = E 2 (g,u 3 (g)); 

e3(g) = E 3 (g,u 3 (g)); e 4 (g) = E 4 (g,u 3 (g)), 

and table II describes the accuracy of the corresponding approximations 

e i(g) = E i(g,ui(g))\ ^(g) = E 2(g,u 3 {g)); e* 3 (g) = El(g,u 3 {g)), 

where E k means the right-hand side of (115). The errors e k (g) = e k {n,g) and £ k (g) = £ k (n,g) depend on 
the value of the coupling parameter g S M + and on the level number n G Z + . The uniform accuracy of an 
approximation may be characterized by the maximal error 

e k = sup sup \e k (n,g)\ 

g eR+ neZ+ 

1G 



or, respectively by 



e* k = sup sup \e* k (n,g)\. 

g eR+ neZ+ 



with the coefficients 



From the table I we have 

ei = 2.0%, e 2 = 0.45%, e 3 = 0.84%, £ 4 = 0.50%, 

and from the table If, 

e[ = 0.40%, e* 2 = 0.37%, £3 = 0.65%. 

As is seen, e k < £& , which means that the corrected approximation (115) improves the accuracy of (114). However, 
the error do not monotonically decrease as the approximation order k increases. This should be related to the 
occurrence of local instabilities in the calculational procedure, which can be detected by the stability analysis. 

To analyse the stability, we, first, need to define the trajectory of the approximation cascade, whose points are 
given by (6). In the considered case we find 

lte(/) = / + /X>fcp(l-y) (116) 

B\\ = A w , 
B21 = ^4.10; B22 = -An, 

B31 = A30, B32 = A31, B33 = A32, 
-841=^30, B42 = A 3 i, B43 = A32, B44 = A33. 

For the quasilocal multipliers (35) we obtain 

M/) = l + E^ P (l + ^)(l-y) P • (117) 

The local multipliers (40) can be found from (117) by means of (41). 

The numerical analysis shows that the condition of quasilocal contraction (47) holds true, 

IM/)I<1 (k = 1,2,3,4), (118) 

and, respectively, the condition of quasilocal stability (56) is valid. This means that all approximations (114), with 
k > 1 are closer to the exact values than the zero approximation. But the condition of local contraction (49) does not 
necessary holds for all / G (1, 00) and n G Z + . Because of this the accuracy of approximations may not improve 
at each step. 

The local Lyapunov exponent (70) for the approximation flow is 

X k (.f)=A k (l+ 2 -^±±^(l-^ k . (119) 

Since / G (1, 00) , the sign of (119) is defined by that of A k , 

sgn Afe (/) = sgm4 fe . 
For A k we have the inequalities 

-4^^il4< 
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-0.002496 <A 2 < 0.003001, 

-0.000525 < ^ 3 < -0.000439 
depending on the quantum number n — 0, 1, 2, . . . . Therefore (119) yields 

Ai (/) < 0.007, A 2 (/) < 0.003, A 3 (/) < 0. (120) 

Thus, at the first two steps the local flow exponents (119) become positive for some energy levels. The positiveness of 
the local Lyapunov exponents signifies the occurrence of local chaos. However, at the third step the motion stabilizes, 

since A3 (/) < for all / G (l,oo) and n G Z + . The same inequalities (120) are valid for the images of Afe (/) 

given by (72), that is, for Afe (3) as functions of g G IR+ and n G Z + . 

VII. CONCLUSION 

Perturbation theory can be made convergent by introducing control functions. The reformulation of perturbation 
theory to the language of dynamical theory makes it possible to control convergence of the approximation sequence by 
checking stability conditions. The sequence of perturbative approximations is bijective to the trajectory of the approx- 
imation cascade. The control functions are defined from quasifixed-point conditions. The terms of the perturbation 
sequence are images of quasifixed points. Each quasifixed point can be corrected by embedding the approximation 
cascade into an approximation flow and considering the motion near the given quasifixed point. 

The stability of calculational procedure is controlled by mapping multipliers and Lyapunov exponents. Several such 
characteristics can be introduced each of them being responsible for controlling particular stability properties. The 
quasilocal multipliers control the quasilocal contraction of a mapping on intervals. The local multipliers control the 
local contraction of a mapping at each step. The quasilocal Lyapunov exponents control the quasilocal stability of the 
approximation cascade on intervals. The local Lyapunov exponents control the local stability of the approximation 
cascade at each point of its trajectory. The classical Lyapunov exponents describe the asymptotic stability of the 
approximation cascade. The local Lyapunov exponents for the approximation flow, enveloping the approximation 
cascade, define the local stability of motion near the corresponding quasifixed points. The stability conditions are 
directly related to the character of convergence of the approximation sequence. Because of the possibility to control 
the convergence by means of control functions, mapping multipliers and Lyapunov exponents, the developed approach 
is called the controlled perturbation theory. 

As an illustration of the method, the calculation of energy levels for a one-dimensional anharmonic oscillator 
is accomplished. Higher-order approximations are considered, as compared to our previous papers. The stability 
analysis showed that, in this case, the approximation cascade is quasilocally stable but not locally stable. In other 
words, the approximation sequence is quasilocally convergent but not locally convergent. This results in local chaos 
at the beginning of the cascade trajectory and, respectively, in a nonmonotonic fluctuation of errors for several first 
approximations. The local chaos can, in principle, be suppressed by introducing a damping parameter diminishing 
the cascade velocity. However, here we did not study the mechanism of such a suppression of chaos. We hope to do 
this in future papers. 
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TABLE I 

The accuracy of the approximations ek(n,g) of the controlled perturbation theory for the energy levels of the 
one-dimensional anharmonic oscillator, as compared to the exact numerical values e(n,g) . 
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n 


e{n,g) 


£i{n,g) 


£i{n,g) 


£a{n,g) 


Si{n,g) 







0.50726 


0.006 


io- 4 


io- 4 


io- 4 




1 


1.5357 


0.009 


0.001 


-io- 4 


-io- 4 




2 


2.5909 


0.001 


0.008 


-0.001 


-io- 4 




3 


3.6711 


-0.008 


0.017 


-0.003 


-0.001 


0.01 


4 


4.7749 


-0.019 


0.027 


-0.006 


-0.002 




5 


5.9010 


-0.030 


-0.002 


io- 4 


-io- 4 




G 


7.0483 


-0.040 


-0.001 


0.001 


io- 4 




7 


8.2158 


-0.051 


-0.001 


0.001 


io- 4 




8 


9.4030 


-0.066 


-0.004 


-0.003 


-0.003 







0.63799 


0.57 


0.10 


0.01 


-0.003 




1 


2.0946 


0.49 


0.15 


-0.01 


-0.004 




2 


3.8448 


-0.06 


0.21 


-0.24 


-0.13 




3 


5.7970 


-0.30 


0.20 


-0.36 


-0.20 


0.3 


4 


7.9100 


-0.40 


0.23 


-0.40 


-0.21 




5 


10.167 


-0.51 


-0.09 


0.01 


-0.01 




6 


12.540 


-0.53 


-0.02 


0.05 


0.03 




7 


15.030 


-0.59 


-0.01 


0.03 


0.01 




8 


17.620 


-0.63 


-0.01 


0.03 


0.01 







0.80377 


1.09 


0.21 


0.02 


-0.01 




1 


2.7379 


0.81 


0.25 


-0.03 


-0.01 




2 


5.1780 


-0.11 


0.25 


-0.40 


-0.22 




3 


7.9400 


-0.40 


0.25 


-0.51 


-0.28 


1 


4 


10.960 


-0.55 


0.24 


-0.58 


-0.32 




5 


14.203 


-0.66 


-0.12 


0.03 


-0.003 




G 


17.630 


-0.69 


-0.05 


0.05 


0.01 




7 


21.240 


-0.77 


-0.06 


0.01 


-0.03 




8 


25.000 


-0.80 


-0.05 


-0.002 


-0.03 







3.9309 


1.97 


0.44 


0.04 


-0.03 




1 


14.059 


1.24 


0.39 


-0.05 


-0.02 




2 


27.550 


-0.24 


0.22 


-0.67 


-0.41 




3 


43.010 


-0.61 


0.18 


-0.79 


-0.48 


200 


4 


60.030 


-0.75 


0.19 


-0.82 


-0.48 




5 


78.400 


-0.85 


-0.18 


0.03 


-0.02 




G 


97.900 


-0.88 


-0.10 


0.03 


-0.02 




7 


118.40 


-0.88 


-0.03 


0.05 


0.01 




8 


139.90 


-0.92 


-0.03 


0.03 


-0.01 







18.137 


2.01 


0.45 


0.04 


-0.03 




1 


64.987 


1.26 


0.40 


-0.05 


-0.02 




2 


127.51 


-0.25 


0.21 


-0.69 


-0.43 




3 


199.20 


-0.64 


0.16 


-0.82 


-0.50 


20000 


4 


278.10 


-0.76 


0.18 


-0.84 


-0.49 




5 


363.20 


-0.84 


-0.16 


0.05 


-0.01 




G 


454.00 


-0.95 


-0.16 


-0.04 


-0.09 




7 


548.90 


-0.91 


-0.05 


0.03 


-0.01 




8 


648.50 


-0.93 


-0.03 


0.03 


-0.01 
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TABLE II 



The accuracy of the corrected approximations 
oscillator, compared to the exact values e(n,g) . 



e* k (n,g) for the energy levels of the one-dimensional anharmonic 



9 


n 


e(n,g) 


£*{n,g) 


el{n,g) 


zl{n,g) 







0.50726 


0.005 


io- 4 


io- 4 




1 


1.5357 


0.007 


0.001 


-io- 4 




2 


2.5909 


0.001 


0.008 


-0.001 




3 


3.6711 


-0.006 


0.017 


-0.003 


0.01 


4 


4.7749 


-0.013 


0.027 


-0.006 




5 


5.9010 


-0.021 


-0.002 


io- 4 




G 


7.0483 


-0.028 


-0.001 


0.001 




7 


8.2158 


-0.035 


-io- 4 


0.001 




8 


9.4030 


-0.046 


-0.004 


-0.003 







0.63799 


0.25 


0.08 


0.01 




1 


2.0946 


0.19 


0.11 


-0.01 




2 


3.8448 


-0.11 


0.11 


-0.22 




3 


5.7970 


-0.20 


0.07 


-0.33 


0.3 


4 


7.9100 


-0.21 


0.06 


-0.36 




5 


10.167 


-0.26 


-0.04 


0.004 




G 


12.540 


-0.23 


0.02 


0.04 




7 


15.030 


-0.26 


0.01 


0.02 




8 


17.620 


-0.27 


0.01 


0.02 







0.80377 


0.28 


0.16 


0.01 




1 


2.7379 


0.19 


0.16 


-0.02 




2 


5.1780 


-0.19 


0.03 


-0.35 




3 


7.9400 


-0.24 


-0.04 


-0.44 


1 


4 


10.960 


-0.26 


-0.08 


-0.49 




5 


14.203 


-0.30 


-0.02 


0.01 




G 


17.630 


-0.28 


0.02 


0.03 




7 


21.240 


-0.32 


-0.02 


-0.02 




8 


25.000 


-0.32 


-0.02 


-0.02 







3.9309 


-0.06 


0.20 


0.001 




1 


14.059 


0.01 


0.13 


-0.03 




2 


27.550 


-0.37 


-0.27 


-0.53 




3 


43.010 


-0.37 


-0.34 


-0.62 


200 


4 


60.030 


-0.33 


-0.35 


-0.64 




5 


78.400 


-0.35 


-0.02 


-0.01 




G 


97.900 


-0.33 


IO" 4 


-0.01 




7 


118.40 


-0.30 


0.04 


0.02 




8 


139.90 


-0.32 


0.01 


-0.01 







18.137 


-0.09 


0.20 


-10" 4 




1 


64.987 


-0.001 


0.13 


-0.03 




2 


127.51 


-0.39 


-0.28 


-0.54 




3 


199.20 


-0.39 


-0.37 


-0.65 


20000 


4 


278.10 


-0.34 


-0.37 


-0.65 




5 


363.20 


-0.33 


-0.01 


0.01 




G 


454.00 


-0.40 


-0.06 


-0.08 




7 


548.90 


-0.32 


0.02 


-0.001 




8 


648.50 


-0.31 


0.01 


-0.003 
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